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Abstract 

We perform a complete "Michel parameter" analysis of all possible helicity structures 
which can appear in the process B — » X c iv^. We take into account the full set of 
operators parametrizing the effective Hamiltonian and include the complete one-loop 
QCD corrections as well as the non-perturbative contributions. The moments of the 
leptonic energy as well as the combined moments of the hadronic energy and hadronic 
invariant mass are calculated including the non-standard contributions. 



1 Introduction 



The experimental and theoretical developments in heavy flavour physics allow us to perform 
a high precision test of the flavour sector. In particular, the enormous amount of data for 
semileptonic B decays in combination with very reliable theoretical methods has opened 
the road for a precision determination of the CKM matrix elements and V u ^,, which are 
known at a relative precision of roughly 2 % and 10 % [1]. 

Aside from testing and extracting its parameters as precisely as possible, a second major 
goal of heavy flavour physics is to look for possible effects beyond the standard model. It 
is genenerally believed that flavour changing neutral currents are a good place to search 
for effects of new physics, since these decays are usually loop-induced and hence sensitive 
to virtual effects from high-mass states. Thus one expects here possibly an effect which is 
sizable compared to the standard-model contribution. 

Semileptonic processes are tree level processes in the standard model and thus the 
relative effects from new-physics contributions are likely to be small. However, a possible 
right-handed admixture to the hadronic current is completely absent in the standard model 
and hence such an effect would be a clear signal for physics beyond the standard model. 

In a recent publication [2] we considered a "Michel parameter analysis" [3] of semilep- 
tonic B decays, where we considered mainly a possible right-handed contribution to the 
hadronic b — > c current. In the present paper we complete the analysis of [2] by extending 
the analysis to all possible two-quark-two-lepton operators. 

There is an extensive literature on a possible non-standard model contributions to 
semileptonic B decays [4, 5, 6, 7, 8]. However, the analysis presented here is different 
in two respects. First of all, our analysis is completely model-independent; however, we 
neglect the lepton masses and hence our analysis would need to be extended straightfor- 
wardly to include e.g. a discussion of a charged Higgs contribution as in [4, 5, 8]. Secondly, 
we consider different observables (i. e. the moments of spectra) which have become avail- 
able only recently through the precise data of the B factories; in this way a much better 
sensitivity to a non-standard contribution is expected. 

In the next section we perform an effective-field-theory analysis of possible new-physics 
contribution, which is kept completely model independent. It turns out that only very 
few operators contribute to the semileptonic b — * c transition. Compared to the usual 
Michel-parameter analysis, well known from muon decay, this effective theory analysis also 
yields order-of-magnitude estimates of the various contributions. 

Based on these effective interactions we recompute the spectra of inclusive semi-leptonic 
b — > c transitons including the new interactons. We make use of the standard heavy 
quark expansion (HQE) and include QCD radiative corrections as well as nonperturbative 
contributions. 

In section 3 we shall perform the HQE including the new-physics operators. In subsec- 
tion 3.1 we compute the QCD radiative corrections for the various helicity combinations 
of the hadronic current. We adopt the kinetic scheme as it has been used for the calcula- 
tion of semileptonic moments in [9] and perform the complete one-loop calculation for the 
new-physics terms. We note that the standard-model calculations have been performed 
already to order a 2 s in [10]. 
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Subsection 3.3 we calculate the nonperturbative contributions of the new physics op- 
erators to order I /ml. Finally, in section 4 we quote our results for the various moments 
which are frequently used in the analysis of semileptonic decays and conclude. 

2 Effective-Field-Theory Analysis of b — ► c£i>£ 

It is well known that any contribution to the effective Lagrangian of some yet unknown 
physics at a high scale A can be written as contributions of operators with dimensions 
larger than four. These operators are SU(3) x SU{2) x 17(1) invariant and suppressed 
by an appropriate power of 1/A. We note that for such an analysis we have to make an 
assumption of the yet not established Higgs Sector: We shall stick with our analysis to the 
single Higgs doublet case; an extension to a type-II two-Higgs doublet as e.g. needed for 
supersymmetry is straightforward. 

The lowest dimension relevant for our analysis is six; the list of relevant operators has 
been given in [11] and we shall use the notations of our previous paper [2]. The quark and 
lepton fields are grouped into 

for the left handed quarks (1) 

for the right handed quarks (2) 

for the left handed leptons (3) 

for the right handed leptons (4) 

where Ql and L L are doublets under SU (2) L and Qr and L R are doublets under an 
(explicitely broken) SU(2) R . Note that we also introduced a right-handed neutrino in 
order to complete the right-handed lepton doubletts. 

The Higgs field and its charge conjugate are written as a 2 x 2 matrix 

tt _ 1 f<Po + iXo V2(j)+ \ ( v 

yft V v/20- 0o - W { ) 

transforming under SU(2)l x SU(2)r. The potential of the Higgs fields leads to a vacuum 
expectation value (vev) for the field 0o- 

The dimension-6 operators fall into two classes, the two-quark operators with gauge 
and Higgs fields and the two-quark-two-lepton operators. In our previous analysis [2] we 
considered only the first class, and the first step towards a full analysis is to also take into 
account the second class. 
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The list of two-quark two-lepton operators with SU(2) L x SU(2) R consists of 1 



0% LL = {Q L Y^ L ){L L Y t L L ) (6) 

vUll = (QLT*r i Q L )(L L T*r i L L ) (7) 

of LfiR = {Q L Y l Q L ){L R Y l L R ) (8) 

{ § RtLL = {Q^iQ R ){L L T i L L ) (9) 
)(0 



OflkflR = (Qr^Qr) (L R Y t L R ) (10) 
while the operators with explicitely boken SU(2) R read 



-PrI.hr = (Qrt^Q^L^Y^) (11) 



Kf L , RR = (QlYiQl) (L R YiT 3 L R ) (12) 

^rr,ll = (Q R Y t r 3 Q R )(L L Y t L L ) (13) 

^rr,rr = (Qr^Qr) {L r Y^L r ) (14) 

Srr,rr = (g R r a r i g R )(L R rV 3 r i L R ) (15) 

Trr,rr = (Q R r a r 3 Y t Q R ) (L R r a r 3 Y t L R ) (16) 

Here we have defined 

Yi (8) Yi = 1 (g) 1, 7/1 ® y, 7^75 <8> 757^, ^ ® (17) 

Note that these operators are not all independent. Furthermore, note that operators with 
the helicity combinations such as (LR)(LR) cannot appear at the level of dimension six 
operators, since additional Higgs fields are required for the helicity flip. 

We shall assume that the right handed neutrino acquires a large majorana mass in 
which case it can be integrated out at some high scale, which we assume to lie well above 
A. In this case SU(2) R is "maximally broken" for the leptons, which means that the 
possible operators always have a projection P_ = (1 — r 3 )/2 and thus only right handed 
interactions involving the right handed charged leptons remain. 

For the case at hand we are interested in the charged current interactions containing a 
b — > c transition. Since we eliminated the right-handed neutrino and helicities are conserved 
for both currents we end up with the conclusion that the charged leptonic current has to 
be left handed. Thus we have only the operators 

o l = {l Ll ,c L ){^ Ll n L ) o 2 = (b RlflCR )(u e , L Yh) (is) 

where £ = e, \i or r, since any helicity changing combination has to originate from dimension- 
8 operators yielding an additional suppression of a factor v 2 /A 2 relative to the dimension-6 
contributions. 



1 In order to have a streamlined notation we suppress all flavour indices in the following 
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However, as has been discussed in our previous paper, helicity violating combinations 
such as (LR) (LL) operators can appear from the two-quark operators with gauge and Higgs 
fields. These operators induce anomalous gauge-boson couplings which are suppressed by 
a factor v 2 /A 2 . They originate from two-quark operators, which are (at the scale of the 
weak bosons) 

Of L = Q L lQ L (19) 
0f L = Q L UQL (20) 

with 

LP = H (iD»H) ] + (iD^H) (21) 
L3 = Hr 3 (iD»H)^ + (iD»H) t 3 H^ (22) 

The terms proportional to r 3 have once again been included to break the custodial sym- 
metry explicitely. 

In the same spirit we define RR-operators 

0^ R = Q R RQ R (23) 

{ r R = Qr{ts,R}Q r (24) 

O l £t = iQ R [r 3 ,M Q R (25) 

{ RR = Q R r 3 Rr 3 Q R (26) 

with 

RP = (iD»H) + {iD^H f H (27) 

Using an odd number of Higgs fields we can construct invariant LR operators. For our 
analysis the relevant operators are 

{ lr = Ql KuB^) HQ r + h.c. (28) 
{ lr = Ql KuWn H Q R + h.c. (29) 
Og> = Q L (iD^H) iD» Q R + h.c. (30) 

After spontaneous symmetry breaking the LL and RR operators contain anomalous 
quark-boson couplings of the order of magnitude v 2 / A 2 . For the LR operators the field 
strenghts of the weak bosons appear, inducing an additional factor of a quark momentum 
p, and hence the order of magnitude is pv/A 2 ~ m q v/A 2 . 

At the scale of the bottom quark theses anomalous coupling terms have the same power 
counting as the two-quark two-lepton operators: Integrating out the weak bosons, their 
propagator together with the gauge couplings reduce to a point like interaction proportional 
to g 2 /M$y = 1/v 2 . Combining this with the order-of-magnitude of the anomalous coupling 
of LL and RR v 2 / A 2 we find that at the scale of the bottom mass these contributions 
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scale in the same way as the two-quark-two-lepton operators directly induced at the high 
scale A. For the case of LR the additional momentum p of the quark is of the order of its 
mass, and hence is as well of the order v 2 /A 2 , possibly further suppressed by a small quark 
Yukawa coupling. 

This conclusion may be altered in a two-higgs doublet model in the case of large tan /?, 
i.e. of a large ratio of Higgs vacuum expectation values. After integrating out the heavy 
degrees of freedom tan/3 will play the role of a coupling constant which then may be 
enhanced by a large value. In e. g. [8] such a scenario has been considered, where a sizable 
value of tan/3 overcomes the supression of the factor (m^mi,) /m 2 H+ in the amplitude; in 
this case also scalar contributions to the leptonic current have to be taken into account. 

To this end, the parametrization introduced in [2] remains valid also in the general case, 
if only dimension-6 operators are included. Thus the effective Hamiltonian reads 

^cff = Jq^Jfi (31) 

where Jj 1 = e7 M P_z/ e is the usual leptonic current and J/^ is the generalized hadronic 
b — > c current which is given by 

Jh,ix = c L cy^P-b + c R c^^P+b + g L ciD^PJb + g R ciD^P+b (32) 
+ d L id u (cia llv P-b) + d R id" ' (cia^P+b) , 

where P± denotes the projector on positive/negative chirality and is the QCD covariant 
derivative. Note that the term proportional to cl contains the standard- model contribution 
as well as a possible new-physics contribution and c R may now also contain a contribu- 
tion from a two-quark-two-lepton operator induced at the high scale A. The gauge part 
ig 3 A®\ a /2 of the QCD covariant derivative gives rise to a new quark-quark-gluon-boson 
vertex. 



3 Operator Product Expansion 

The operator product expansion (OPE) for inclusive decays has become textbook material 
[12]. For the case of inclusive semileptonic decays the OPE is formulated for the T-product 
of the two hadronic currents 

T, u = J rf 4 xe-^ m ^(5(p)|^(x)r M c( a :)c(0)r^(0)| J B(p)) (33) 

where V is the combination of Dirac matrices and derivatives given in (32), v = p/M R 
is the four velocity of the decaying B meson and q is the momentum transferred to the 
leptons. The quantity is expanded in inverse powers of the scale of the order mj, where 
mj is the the heavy quark mass. Technically this procedure is an OPE for the product of 
the two currents. 
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Figure 1: Real Corrections 
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Figure 2: Virtual Corrections 



The standard- model calculation has been performed at tree level up to order 1/mf and 
it turns out that the non-perturbative corrections are small. The radiative corrections have 
been computed to order a s , (3oa 2 s and recently also to order a 2 s for the leading (i.e. the 
parton model) term [10] and to order a s for the term of order l/m 2 involving /i 2 . 

In the following we shall consider the perturbative and non-perturbative contributions 
to the OPE, performed with the modified current (32). We shall compute the complete 
one-loop contributions as well as the leading non-perturbative corrections proportional to 
/i 2 and /Iq. 

3.1 QCD Corrections and Renormalization Group Analysis 

The calculation of the QCD radiative corrections has been performed in [13] and the results 
in the kinetic scheme have been given in [9] for the semileptonic moments in the standard 
model. In order to perform an analysis of possible non-standard contributions we have to 
calculate the QCD radiative corrections for the current (32) to order a s . Thus we have to 
evaluate the Feynman diagrams shown in fig. 1 and 2 for the real and virtual corrections 
respectively. Note that the scalar current (i.e. the terms proportional to Ql/r) induces 
new vertices shown in the Feynman diagrams at right. The real and virtual corrections are 
individually IR-divergent. We regulated the IR-divergence by introducing a gluon mass 
which drops out upon summation of the real and virtual correction being IR-convergent. 
In the calculation of the virtual corrections the wave function renormalizations of the b 



6 



and c quark field also have to be included. 

The total amplitude consists of the sum of the standard- mo del contribution and the 
one from the new-physics operators. Since the new-physics piece is of order 1/A 2 , we 
shall include only the interference term of the standard model with this contribution. The 
square of the new-physics term is already of order 1/A 4 and has to be neglected, since we 
compute only up to this order. Thus we compute 

L ■((clv\c L H™\b)(clv\H cS \by + {dv\?iri\b){dv\c L H^\b)*)d<l> PS (34) 



dr 



2m h 



where d0ps is the corresponding phase space element and 

■(c^yP-fe) (e 7 M P-^ e ) 



n/SM 



4G F V cb , 



V2 



is the standard-model effective Hamiltonian, which has the same helicity structure as the 
new-physics contribution proportional to Cl- 

The relevant Feynman rules for the new-physics operators at tree level can be read off 
from (32); note that the terms involving g L and g R yield a boson-gluon-quark-antiquark 
vertex in order to maintain QCD gauge invariance. 

It is well known that the left- and right-handed currents do not have anomalous dimen- 
sions and hence the parts of (32) with cl and cr are not renormalized. However, the scalar 
and tensor contributions have anomalous dimensions and hence we need to normalize these 
operators at some scale and run them down to the scale of the bottom quark. 

To this end, we have to calculate the anomalous-dimension matrix of these currents to 
set up the renormalization group equation. It can be obtained from the requirement that 
the physical matrix elements must not depend on the renormalization scale //: 

V = ^T—(^e\H cfi \b} (35) 
am fi 

Inserting the OPE for the Hamiltonian we get: 



{du e \n cS \b) 



with 



C 



V2 

4C F Kb 
V2 

( 9 l\ 

9r 
d L 
d R 



(ciu e \ [c L {c llx P-b){e^P-V e ) + c R (c llx P + b)(e^P-V e )} \b) 



C ■ (c£v e \0\b), 



(36) 



c 



L 

m b 
R 



o 



VrJ 



( {ci%P_b){e^P_v e ) \ 
{ci%P + b){e^P_v e ) 
(id»(cia^P-b))(e^P_v e ) 

(id»(cia^P + b))(erP-V e ) 
(m b c7 M P_6)(e7 M P_w e ) 
(m fe c7 M P + 6)(e7^P_v e ) 
(m c c7 M P_6)(e 7 ^P_t; e ) 
V (m c c 7M P + 6)(e 7 ^P_t; e ) / 



(37) 
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where the operators O are of dimension seven. 

In the following we consider the renormalization group mixing of these dimension-four 
operators. The calculation of the one-loop anomalous dimension is standard. We define 
the anomalous dimension matrix 7 by: 



dC 
din// 



= 7 T (/i) C 



(38) 



and compute 7 from the divergencies of the renormalization constants in the usual way. 
We obtain 



7 T M 



3tt 



/o 




















o\ 
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V 



(39) 



The renormalization group equation for the Wilson coefficient is 



v (91n/x ' ^ v s da s 
We seek a solution of this equation with the initial conditions 

Cl/rW = = c™J R (A) , 



L/R\ 



(40) 



(41) 



since the matching of the left- and right-handed currents is performed by fixing the coef- 
ficients Cl and Cr and all additional contributions are only due to renormalization group 
running. Inserting the one-loop expressions we obtain 

cl/r{h) = CL/i?(A) 
9l/r(v) = 9l/r(A) 



d L /R{lA 



{g L/R (A) + d L/R (A)) ( ^) - g L/R (A) 



"L/R 




(42) 



^L/R 



One may reexpand (42) using the one-loop expression for the strong coupling constant and 
obtain the logarithmic terms of the one-loop calculation. However, the straight-forward 
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one-loop calculation also yields finite terms, which depend on the choice of the renormal- 
ization scale fi. It is well known that in order to fix this dependence on the renormalization 
scale, one would need to include the running at two loops, which, however, goes beyond 
the scope of the present paper. Rather we shall fix this scale to be /i = mj,, which is the 
relevant scale of the decay process, assuming that the full NLO calculation would fix a 
scale of this order. 

The advantage of this procedure is that the kinematic effects, which lead to a distortion 
of the spectra and thus have an impact on the moments, are given by these finite terms. 
We expect that a full NLO calculation will lead to very similar results. 

3.2 Mass Scheme 

The calculation of the process is usually set up with pole masses of the particles, which is 
known not to be a well defined mass. The problems manifest themselves by abnormally 
large radiative corrections when the pole mass scheme is used. It has been discussed 
extensively in the literature that an appropriately defined short- distance mass is better 
suited for the OPE calculation of an inclusive semileptonic rate. 

In the present analysis we will use the kinetic mass scheme, where the mass is defined 
by a non-relativistic sum rule for the kinetic energy [14] . At one-loop level the kinetic mass 
is related to the pole mass by 



where \if is a factorization scale for removing contributions below from the mass-definition. 
The factorization scale is set to 1 GeV since this is the typical energy release in the process. 
This low renormalization scale is in fact the reason why the MS scheme is inappropriate. 

The ratio p = m 2 c /m^ is rather stable under the choice of schemes (provided that the 
same scheme is chosen for both and m c ) and thus the choice of the mass scheme enters 
only through the dependence of the rates. It is well known from the calculation in the 
standard model that the 0(a s ) corrections from the relation of the kinetic mass with the 
pole mass 



compensate to a large extent the radiative corrections to the rates computed in the pole 
scheme, leaving only small QCD radiative corrections. It turns out that this also is the 
case in our calculation including the anomalous couplings. 

3.3 Non-perturbative Corrections 

The nonpertubative corrections at tree level, including the modified current (32), have 
been studied in [15], however, these results have not yet been published, and hence we 
shall quote these results in the following. 




(43) 




(44) 
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The calculation of the nonperturbative corrections at tree level requires to compute 
the Feynman diagram shown in fig. 3, where the one-gluon graph is needed to obtain the 
matching coefficient of the chromo-magnetic moment operator. 

The calculation is standard and yields somewhat lengthy results, thus we defer the 
presentation of these expressions to the appendix. 



w- w- w~ 





b b b g 

Figure 3: Nonpertubative corrections at tree level 



To the order we calculate, the nonperturbative effects are parametrized by the kinetic 
energy [i 2 and the chromomagnetic moment /Iq, which are small quantities compared to 
the 6-quark mass. Inserting the values extracted from b —> c semileptonic decays we find 

2 2 

/,;: /,;/ - 0.02. (45) 



ml 3 m 



Hence the non-perturbative corrections are tiny compared to the leading terms, as long 
as there are no abnormally large coefficients or the leading term vanishes. As it has been 
investigated in [15] this as well holds true for the new-physics contributions parametrized 
by (32). 



4 Results and Discussion 

We have evaluated the new-physics contributions to the various moments of the leptonic 
and hadronic energy and the hadronic invariant mass spectra. We have included tree-level 
partonic and l/mf corrections as well as the QCD radiative corrections at one- loop with 
a renormalization group treatment as described in the last section. The hadronic energy 
and the hadronic invariant mass of the decay products can be written as 

E Ha d = v ■ (p B - q) = m B - v ■ q 
SHad = (pb - qf = m 2 B - 2m B vq + q 2 , 

where m# and ps = tubV are the mass and the momentum of the B meson and q is the 
momentum of the leptonic system. The S-meson mass can be expanded as 

2 2 

m B = m b + A + ^-t^ + ... . (47) 
2m h 
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Thus it is possible to relate the hadronic variables in (46) to the partonic ones 

^ , (48) 

s = — o = 2 = l-2v ■ q + q , 

m b m b 

where is the the b quark momentum. In the following we shall quote the results in terms 
of the partonic variables (48). 

We have also included a cut on the charged lepton energy since such a cut has to be 
used in the experimental analysis. The results may be obtained as FORTRAN code from 
the authors. In order to have a qualitative discussion of the results, we give the results 
in tables 1-3 for various moments without an energy cut for the charged lepton energy E\ 
and with a cut of 1 GeV for this quantity. In table 1 we list the results for the moments 

L n = ± [ dE, K ^~ , (49) 
1 o JE cut dEi 

in table 4 we consider the scale dependence of the L n , and in the tables 2 and 3 we quote 

If , A f , A ,„ „„■ d 3 r 



= ±- I dE t I ds dE (s - pf 

J- JEc U t J 



o JE cut J dE ds dEi 

with p = ml /ml, where the normalization 

To = ° l ^f b (1 - 8P - 12p 2 lnp + 8p 3 - p 4 ) (50) 

is given in terms of the partonic rate. Note that we perform the calculations in the kinetic 
scheme, and we also insert the value of the kinetic mass in the normalization IV 

The entries in the tables contain the coefficients corresponding to the expansion of the 
various moments: 

L n = c\L^ + c L c R L^ + c L d L L^) (51) 
+ c L d R L^n) + CL g L L^ + c L g R L^ 

H l3 = c\H^ + c L c R H^ CR) + c L d L H^ (52) 
+ c L d R H^ dR) + c L g L H^ + c L9R H^ 

where all the coefficients have an expansion in a s and in l/m& 

2 2 

r(ciC2) _ T ( c l c 2; m °. a °) i f^ir r(cic 2 ;ml,a° s ) T (cic 2 ; ,a°) a s r (cic 2 ; mj,^) , 
-k„ — Lin ~\ ^Li n -+- 7;Li n -+- ■ ■ ■ -\ Li n -+- ■ ■ ■ 

ml 3 ml n 

2 2 

jj(cic 2 ) _ jj(cic 2 ;m°,a s ) /V^-(cic 2 ; mg,a°) /^g jj(ciC2; mg,a°) _ _ _ «s^-(cic 2 ; mg,aj) _ _ _ 

ij ~ ij ///•; 3 ///•; ij 7T ij ' 
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where we have only shown the terms which we have calculated. The values of couplings 
to be inserted in (51) and (52) are the ones at the high scale A. Furthermore, for the 
numerical analysis we use m^ m (l GeV) = 4.6 GeV and p = m 2 c /ml = 0.0625. 

The results of the calculations are displayed in tables Tab. 1-5 in the appendix. Table 
1 contains the results for the leptonic moments normalized to the total leptonic rate at 
tree level (50) for all lepton energies and for a cut of 1 GeV on the lepton energy. It turns 
out that the radiative corrections to the scalar and tensor admixtures are sizable, i. e. the 
a s /ir coefficients are large. In addition, these coefficients have the opposite sign as the tree 
level piece, and hence a substantial reduction of the tree result is expected. 

Table 4 contains the sum of the tree level and the a s contributions using the one-loop 
expression for the running coupling a s . As discussed above, the full NLO expressions for 
the scalar and tensor couplings are not available yet and hence a residual scale dependence 
remains. We expect the scale to be of the order of m b and hence we evaluate the expressions 
for /i = m b /2, m b and 2m b . For c| as well as for c^cr the scale dependence is weak and 
originates from yet unknown NNLO effects. Due to the large a s /n coefficients the scale 
dependence for the tensor couplings is sizable, while it is huge for the scalar couplings, 
since the tree contribution is almost cancelled by the radiative correction. A full NLO 
calculation will very likely not improve this situation and hence we have to conclude that 
we will not have a good sensitivity to the tensor couplings and practically no sensitivity to 
the scalar couplings, at least for the lepton-energy moments. 

The coefficients of the nonperturbative contributions at tree level are in general of 
similar size as the ones of the a s corrections. Since a s /ir ~ p^/ml, the non-perturbative 
corrections are of similar importance as the radiative ones. However, the leptonic moments 
are all dominated by the tree-level contribution and hence the radiative as well as the 
nonperturative corrections to the moments are small. 

Tables 2 and 3 contain the various hadronic moments computed without and with a 
cut on the lepton energy. For the % — moments we have to draw the same conclusion 
as for the leptonic moments: The scalar and tensor couplings have large and opposite- 
sign coefficients compared to the tree level piece; this leads in the same way to a sizable 
reduction of the tree level result as well as to a large scale dependence, which is shown in 
table 5, where the result up to order a s is shown. 

Clearly the moments with i > do not have a tree level contribution at the partonic 
level since the tree-level partonic rate is proportional to the mass shell delta function 
5(s — p). For these moments the leading contributions are at order a s or l/m£. Hence 
their dependence on the scale is given by the dependence of a s . However, here the radia- 
tive corrections are small compared to the non-perturbative ones. The non-perturbative 
corrections at tree level contain also derivatives of the mass shell delta function 5(s — p), 
where at leading order l/m b the maximum number of derivatives is two. Due to this, the 
first and second % moments are of order l/mf; higher moments with % > 2 will only have 
contributions of order 1/mf or higher. 

The sensitivity to a possible new-physics contribution is mainly limited by the precision 
of the standard-model calculation. Current analyses use up to the second moments in both 
the leptonic energy and the invariant mass squared. The highest moments included in the 
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standard-model analyses are (roughly) sensitive to terms of the order l/mf which is the 
highest order in the l/m& expansion included in the fit. The size of this terms together 
with the size of the a 2 corrections may serve as a conservative estimate of the uncertainties 
of the standard model calculation, which at the end determines the sensitivity to a possible 
new-physics contribution. Furthermore, an inclusion of higher moments in the fit, including 
the new contributions, (in particular with i > 2) needs the calculation of the l/mf terms 
for the new-physics contributions. As the impact of such hadronic mass moments to the fit 
is small we did not include a table of them in this paper, but the results of the calculation 
can be obtained by the authors in a Fortran or Mathematica file. 

5 Summary and Conclusions 

This work completes the analysis of possible new physics effects in inclusive semileptonic 
B decays. Starting from a general ansatz for anomalous couplings in semileptonic decays 
we compute the effects on leptonic and hadronic moments which are used in the analysis 
of inclusive semileptonic decays. 

As far as the leptonic moments are concerned, the QCD radiative corrections turn out to 
be as important as the nonperturbative ones. We have presented the complete expressions 
to order a s and to order l/m 2 including the new-physics pieces. 

This holds also true for the hadronic energy moments. However, the hadronic mass 
moments (taken with respect to m 2 ) do not have a tree-level contribution. Hence the 
nonperturbative corrections of order 1/m 2 as well as the terms of order a s are the leading 
contributions in the heavy quark expansion. It turns out that, numerically, the non- 
perturbative contributions are in general dominant. 

For the leptonic moments and for the % — hadronic moments the radiative correction 
for the scalar and tensor couplings turn out to be sizable. This leads to a substantial 
reduction of the moments and in the case of the scalar coupling to an almost cancellation 
between tree level and the radiative corrections, which induces a large scale dependence. 
Hence the sensitivity to scalar and tensor couplings of the moments is limited. However, 
the moments with i ^ appear first at order a s and have a resonable sensitivity to scalar 
and tensor couplings. 

The standard analysis in semileptonic decays is to perform a combined fit of V c ^, the 
quark masses and the HQE parameters, usually up to order 1/mf. We propose to use the 
results given here to include the anomalous couplings induced by possible new physics into 
such a fit. 

The effective-theory analysis indicates that a right-handed admixture could be the 
largest effect. Since the radiative corrections to the right-handed currents are completely 
known to NLO and the size of the coefficients indicate a good sensitivity to the anomalous 
coupling cr this coefficient should be the first to be searched for in a moment analysis. 
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A Nonperturbative corrections to the non-standard 
currents 

In this appendix we show the results for the new-physics contributions based on (32). We 
have calculated the interference term of the standard-model contribution with the current 
(32) and list the terms proportional to the coupling constants. We obtain 
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-12(y-iy-24py 2 m b 



2(4 2 / 2 -9y + 3)pV lo 2 , fi 2 \ pi 
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B Tables 



n 


c| c L c R c^gf/j c^di Cidjj 


a 

3 


1.0000 -0.6685 0.2212 0.5400 0.3315 -0.6597 
0.3072 -0.2092 0.0613 0.1372 0.0977 -0.2307 
0.1030 -0.0708 0.0188 0.0388 0.0314 -0.0845 
0.0365 -0.0252 0.0062 0.0118 0.0107 -0.0319 


13 

CP Er « 1 

"3 $^82 

^ 3 


-0.5000 0.3342 -0.0017 0.1703 -0.1652 0.3288 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0858 -0.0590 0.0365 0.1146 0.0261 -0.0702 
0.0730 -0.0503 0.0210 0.0575 0.0214 -0.0637 


§ 

c cK 1 

^ cu 1 

"J* 8 2 
3 


-1.9449 4.9934 1.0232 1.5624 -2.1536 3.7106 
-0.9625 1.8578 0.3253 0.6011 -0.7986 1.5873 
-0.4495 0.7237 0.1124 0.2427 -0.3081 0.6840 
-0.2052 0.2902 0.0410 0.1008 -0.1220 0.2966 


•S 
8 1 

\ 2 
e 3 


0.3125 0.8009 -2.6592 -8.8212 -2.1497 4.3637 
0.0908 0.2284 -0.7171 -2.3141 -0.5594 1.4880 
0.0276 0.0739 -0.2174 -0.6843 -0.1660 0.5394 
0.0085 0.0260 -0.0711 -0.2189 -0.0538 0.2039 


« 
1 

fH 

£ 2 
3 


0.8148 -0.5617 0.1621 0.3586 0.2631 -0.6161 
0.2776 -0.1919 0.0520 0.1089 0.0867 -0.2232 
0.0979 -0.0678 0.0172 0.0340 0.0296 -0.0831 
0.0356 -0.0246 0.0059 0.0109 0.0104 -0.0317 


13 

> ^^82 

OH 3 
1 — 1 


-0.4504 0.3225 0.0433 0.3440 -0.1479 0.3631 
0.0087 -0.0021 0.0564 0.2247 0.0031 0.0059 
0.0874 -0.0594 0.0377 0.1194 0.0267 -0.0691 
0.0733 -0.0504 0.0213 0.0583 0.0215 -0.0635 


A 
&q j=_ 1 
«^ 2 
3 


-2.1029 4.6903 0.8592 1.4595 -2.0451 3.7102 
-0.9883 1.8078 0.2989 0.5845 -0.7805 1.5871 
-0.4540 0.7149 0.1078 0.2398 -0.3049 0.6840 
-0.2060 0.2886 0.0401 0.1003 -0.1214 0.2966 


•S 
8 1 

^ 2 
e 3 


0.2640 0.5740 -1.8506 -5.9374 -1.3992 3.9213 
0.0828 0.1930 -0.5920 -1.8692 -0.4440 1.4126 
0.0262 0.0679 -0.1964 -0.6098 -0.1467 0.5260 
0.0083 0.0249 -0.0674 -0.2058 -0.0504 0.2014 



Table 1: Tree level and a s /rc coefficients of the leptonic moments without Ei cuts and with 
a cut Ei>l GeV. Note that we have redefined cIl/r = me cIl/r and gL/R = tub gi/R with 
rriB = 5.279 GeV in order to tabulate dimensionless quantities. 
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i 


j 


c 2 


c L c R 


c L 9l 


Cl9r 


cl<Il 


CLd R 













1.0000 


-0.6685 


0.2212 


0.5400 


0.3315 


-0.6597 




a 

n 





1 


r\ a c~\c~\ r\ 

0.4220 


-0.2500 


0.0961 


0.2556 


0.1217 


-0.2559 




parti 





2 


0.1832 


-0.0964 


0.0429 


0.1219 


0.0461 


-0.1021 







3 


0.0815 


-0.0383 


0.0196 


0.0586 


0.0180 


-0.0418 






i > 


j 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 












-0.5000 


0.3342 


-0.0017 


0.1703 


-0.1652 


0.3288 









1 


-0.5000 


0.3342 


-0.100 


-0.2229 


-0.1652 


0.3288 




CP 





2 


-0.2902 


0.1836 


-0.0773 


-0.2119 


-0.0899 


0.1840 


i — i 





3 


-0.1382 


0.0837 


-0.0448 


-0.1348 


-0.0406 


0.0853 


i 


O 


1 





-0.5780 


0.4185 


-0.2038 


-0.5937 


-0.2091 


0.4025 


Tree 




1 


1 


-0.1584 


0.1172 


-0.0695 


-0.2158 


-0.0585 


0.1129 


CM S= 

3. 


1 


2 


-0.0283 


0.0280 


-0.0217 


-0.0718 


-0.0143 


0.0258 






2 





0.1609 


-0.0728 


c\ Aon/ 1 

0.0386 


0.1159 


0.0337 


-0.0809 






2 


1 


0.0735 


-0.0302 


0.0180 


0.0561 


0.0138 


-0.0343 








n 


n nnnn 


n nnnn 


n nnnn 


n nnnn 


n nnnn 


— n nnnn 












— 1.9449 


a r\ r\ o a 

4.9934 


1.0232 


1.5624 


-2.1536 


3.7106 




cK 

o 





l 


-0.3850 


1.2777 


c\ a r\r\ n 

0.4097 


0.4782 


— n.5223 


0.9700 







2 


-0.0302 


0.2833 


0.1576 


0.1391 


—0.1109 


0.2254 




o 
o 





3 


0.0298 


0.0342 


0.0578 


0.0350 


-0.0146 


0.0347 






1 





0.3143 


-0.6395 


-0.1100 


-0.2167 


0.2027 


-0.4360 




a. 


1 


1 


0.1195 


-0.2561 


-0.0529 


-0.0925 


r\ r\ p—t a A 

0.0744 


-0.1709 






1 


2 


0.0466 


-0.1059 


-0.0254 


-0.0405 


0.0282 


-0.0689 






i > 1 


J 


0.0000 


r\ r\ r\r\ r\ 

0.0000 


r\ r\ r\r\ r\ 

0.0000 


r\ r\r\ r\r\ 

0.0000 


0.0000 


r\ r\ r\ r\r\ 

0.0000 












0.3128 


0.8007 


-2.6592 


-8.8212 


-2.1497 


4.3637 









1 


0.1631 


0.3441 


-1.2391 


-4.1901 


-0.8839 


1.8575 









2 


0.0910 


0.1477 


-0.5850 


-2.0067 


-0.3694 


0.8017 


yd 

CD 







3 


0.0526 


0.0632 


-0.2793 


-0.9681 


-0.1568 


0.3505 


o 

CJ 




1 





0.0901 


-0.0363 


0.0028 


0.0176 


0.0032 


-0.0095 


l> 




1 


1 


0.0470 


-0.0178 


0.0014 


0.0093 


0.0015 


-0.0046 


VI 

a 




1 


2 


0.0251 


-0.0090 


0.0007 


0.0050 


0.0007 


-0.0023 






2 





0.0091 


-0.0033 


0.0001 


0.0015 


0.0002 


-0.0008 






2 


1 


0.0053 


-0.0019 


0.0000 


0.0009 


0.0001 


-0.0004 






3 





0.0018 


-0.0006 


0.0000 


0.0003 


0.0000 


-0.0001 



Table 2: Tree level and a s /ir coefficients of the hadronic moments without E\ cuts. The 
partonic tree- level moments for i > 1 are all zero. Note that we have redefined dL/R = 
iTiBdi/R and gL/R = rnB9hiR with uib = 5.279 GeV in order to tabulate dimensionless 
quantities. 



17 



i 


j 


c 2 


c L c R 


c L 9l 


Cl9r 


cl<Il 


CLd R 













0.8148 


-0.5617 


0.1621 


0.3586 


0.2631 


-0.6161 




a 

n 





1 


0.3341 


-0.2037 


0.0682 


0.1676 


0.0922 


-0.2365 




parti 





2 


0.1411 


-0.0761 


0.0295 


0.0789 


0.0332 


-0.0933 







3 


0.0612 


-0.0293 


0.0131 


0.0375 


0.0123 


-0.0378 






i > 


j 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 












-0.4504 


0.3225 


0.0433 


0.3440 


-0.1479 


0.3631 









1 


-0.4505 


0.2921 


-0.0597 


-0.0843 


-0.1329 


0.3332 




CP 





2 


-0.2673 


0.1561 


-0.0532 


-0.1300 


-0.0695 


0.1841 


i — i 





3 


-0.1337 


0.0706 


-0.0327 


-0.0935 


-0.0308 


0.0859 


i 


O 


1 





-0.5424 


0.3590 


-0.1687 


-0.4845 


-0.1685 


0.3887 


Tree 




1 


1 


-0.1639 


0.1022 


-0.0598 


-0.1852 


-0.0478 


0.1115 


CM S= 

3. 


1 


2 


-0.0417 


0.0262 


-0.0204 


-0.0678 


-0.0126 


0.0273 






2 





0.1203 


-0.0547 


0.0258 


0.0742 


0.0223 


-0.0729 






2 


1 


0.0538 


-0.0221 


0.0118 


0.0355 


0.0087 


-0.0306 








n 


n nnnn 


n nnnn 


n nnnn 


n nnnn 


n nnnn 


n nnnn 












—2.1029 


4.6903 


0.8592 


1.4595 


—2.0451 


3.7102 




cK 

o 





l 


-0.4609 


1.2205 


0.3461 


0.4476 


-0.5005 


0.9855 







2 


-0.0660 


0.2921 


0.1348 


0.1332 


—0.1119 


0.2391 




o 
o 





3 


0.0131 


0.0538 


0.0507 


0.0363 


-0.0194 


0.0439 






1 





0.3074 


-0.5095 


-0.0803 


-0.1804 


0.1654 


-0.4093 




a. 


1 


1 


0.1171 


-0.1971 


-0.0381 


-0.0751 


0.0583 


-0.1590 






1 


2 


0.0458 


-0.0789 


-0.0180 


-0.0321 


0.0211 


-0.0635 






i > 1 


J 


0.0000 


r\ r\ r\r\ r\ 

0.0000 


r\ r\ r\r\ r\ 

0.0000 


r\ r\r\ r\r\ 

0.0000 


r\ r\r\ r\r\ 

0.0000 


0.0000 












0.2642 


0.5739 


-1.8506 


-5.9373 


-1.3992 


3.9213 









1 


0.1216 


0.2462 


-0.8449 


-2.7806 


-0.5529 


1.6572 









2 


0.0608 


0.1057 


-0.3919 


-1.3149 


-0.2221 


0.7103 


yd 

CD 







3 


0.0323 


0.0455 


-0.1842 


-0.6272 


-0.0907 


0.3086 


o 

CJ 




1 





0.0576 


-0.0231 


0.0018 


0.0101 


0.0018 


-0.0079 


l> 




1 


1 


0.0288 


-0.0108 


0.0009 


0.0052 


0.0008 


-0.0038 


VI 

a 




1 


2 


0.0147 


-0.0052 


0.0004 


0.0027 


0.0004 


-0.0018 






2 





0.0046 


-0.0016 


0.0001 


0.0007 


0.0001 


-0.0006 






2 


1 


0.0026 


-0.0009 


0.0000 


0.0004 


0.0000 


-0.0003 






3 





0.0007 


-0.0002 


0.0000 


0.0001 


0.0000 


-0.0001 



Table 3: Tree level and a s /n coefficients of the hadronic moments with a cut Ei>l GeV. 
The partonic tree-level moments for i > 1 are all zero. Note that we have redefined d L / R = 
iTiBdi/R and gL/R = rnB9hiR with uib = 5.279 GeV in order to tabulate dimensionless 
quantities. 
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[i n 


<4 c L c R Cigfx, c^str c^di Cid R 


> 

o i 

oo 2 
N 3 


1.0253 -0.6037 0.0042 -0.1916 0.1533 -0.2983 
0.3145 -0.1907 0.0028 -0.0552 0.0512 -0.1074 
0.1052 -0.0648 0.0011 -0.0182 0.0176 -0.0397 
0.0372 -0.0231 0.0004 -0.0065 0.0062 -0.0150 


> 

^D 2 

^ 3 


1.0208 -0.6151 0.0441 -0.0474 0.1883 -0.3692 
0.3132 -0.1940 0.0135 -0.0169 0.0604 -0.1317 
0.1048 -0.0658 0.0043 -0.0068 0.0204 -0.0485 
0.0371 -0.0234 0.0015 -0.0028 0.0071 -0.0184 


> 
6 1 

<N 2 

^ 3 


1.0177 -0.6231 0.0715 0.0752 0.2146 -0.4223 
0.3123 -0.1963 0.0208 0.0164 0.0674 -0.1499 
0.1046 -0.0666 0.0065 0.0034 0.0225 -0.0552 
0.0370 -0.0237 0.0022 0.0005 0.0078 -0.0209 



Table 4: Summed up tree level and a s /n coefficients of the leptonic moments without E\ 
cuts for n = 2.3, 4.6 and 9.2 GeV. 



M i J 


c 2 L c L c R c L g L c L g R c L d L c L d R 


> 

o 1 
« 2 

04 3 


1.0253 -0.6037 0.0042 -0.1916 0.1533 -0.2983 
0.4352 -0.2222 -0.0051 -0.0914 0.0486 -0.1024 
0.1906 -0.0845 -0.0049 -0.0440 0.0156 -0.0360 
0.0857 -0.0331 -0.0033 -0.0214 0.0050 -0.0129 


> 

6 1 

^02 
^ 3 


1.0208 -0.6151 0.0441 -0.0474 0.1883 -0.3692 
0.4329 -0.2271 0.0136 -0.0234 0.0628 -0.1322 
0.1892 -0.0866 0.0040 -0.0117 0.0215 -0.0487 
0.0850 -0.0340 0.0010 -0.0059 0.0075 -0.0185 


> 

6 1 

^02 
^ 3 


1.0177 -0.6231 0.0715 0.0752 0.2146 -0.4223 
0.4312 -0.2305 0.0269 0.0335 0.0734 -0.1545 
0.1883 -0.0880 0.0104 0.0151 0.0258 -0.0582 
0.0845 -0.0347 0.0041 0.0069 0.0093 -0.0226 



Table 5: Summed up tree level and a s /n coefficients of the non- zero-tree-level hadronic 
moments without E t cuts for \x = 2.3,4.6 and 9.2 GeV. 



19 



References 

[1] See e. g. M. Neubert, Contribution to the 23rd International Symposium on Lepton- 
Photon Interactions at High Energy (LP07), Daegu, Korea, 13-18 Aug 2007 
[arXiv:hep-ph/0801.0675]. 

[2] B. M. Dassinger, R. Feger and T. Mannel, Testing the Left- Handedness of the b — > c 
Transition, Phys. Rev. D 75, 095007 (2007) [arXiv:hep-ph/0701054]. 

[3] L. Michel, Interaction Between Four Half-Spin Particles and the Decay of the ix- Meson, 
Proc. Phys. Soc. A 63, 514 (1950). 

[4] Y. Grossman and Z. Ligeti, The Inclusive B — > tvX Decay in Two Higgs Doublet 
Models, Phys. Lett.B332, 373 (1994) [arXiv:hep-ph/9403376]. 

[5] Y. Grossman and Z. Ligeti, Transverse Tau Polarization in inclusive B — > tvX decays, 
Phys. Lett.B347, 399 (1995) [arXiv:hep-ph/9409418]. 

[6] M. B. Voloshin, Bound on V + A Admixture in the b — > c Current from Inclusive 
vs. Exclusive Semileptonic Decays of B Mesons, Mod. Phys. Lett. KYI, 1823 (1997) 
[arXiv:hep-ph/9704278]. 

[7] T. G. Rizzo, Right-Handed Currents in B Decays Revisited, Phys. Rev.D58, 055009 
(1998) [arXiv:hep-ph/9803385]. 

[8] U. Nierste, S. Trine and S. Westhoff, Charged-Higgs effects in a new B — > D tau 
nu differential decay distribution, Phys. Rev. D 78, 015006 (2008) [arXiv:0801.4938 
[hep-ph]]. 

[9] P. Gambino and N. Uraltsev, Moments of Semileptonic B Decay Distributions in the 
l/m b Expansion Eur. Phys. J. C 34, 181 (2004) [arXiv:hep-ph/0401063]. 

[10] A. Czarnecki and K. Melnikov, Semileptonic b — > u decays: Lepton invariant mass 
spectrum, Phys. Rev. Lett. 88, 131801 (2002) [arXiv:hep-ph/0112264]. 

[11] W. Buchmiiller and D. Wyler, Effective Lagrangian Analysis of New Interactions and 
Flavour Conservation, Nucl. Phys. B 268, 621 (1986). 

[12] A. V. Manohar and M. B. Wise, Heavy Quark Physics, Camb. Monogr. Part. Phys. 
Nucl. Phys. Cosmol. 10 (2000) I. 

[13] M. Jezabek and J. H. Kuhn, Lepton Spectra From Heavy Quark Decay, Nucl. Phys. 
B 320, 20 (1989). 

[14] I. I. Y. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Sum Rules for 
Heavy Flavor Transitions in the SV Limit, Phys. Rev. D 52, 196 (1995) [arXiv:hep- 
ph/9405410]. 



20 



[15] B. M. Dassinger, Semileptonische B-Zerfdlle als Test des Standardmodells, Diploma 
Thesis, University of Siegen, (2006) [in German] 

R. Feger, Test der V-A-Struktur von Quark- Strdmen im inklusiven Zerfall B — > 
X c e~u e , Diploma Thesis, University of Siegen, (2005) [in German] 



21 



